Regular and chaotic vortex core reversal by a resonant perpendicular field 
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Under the action of an alternating perpendicular magnetic field the polarity of the vortex state 
nanodisk can be efficiently switched. We predict the regular and chaotic dynamics of the vor- 
tex polarity and propose a simple analytical description in terms of a reduced vortex core model. 
Conditions for the controllable polarity switching are analyzed. 
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I. INTRODUCTION 

Investigation of the magnetization dynamics at the 
nanoscale is a key task of the modern nanomagnetism.^ 
One of the typical topologically nontrivial magnetic con- 
figurations of a nanoscaled magnet is a magnetic vor- 
tex, which can form a ground state configuration of 
submicron-sized magnetic disk-shaped particles (nan- 
odots). Such a vortex is characterized by a curling 
divergent -free in-plane configuration with magnetization 
tangential to the edge surface of the nanoparticle.^ The 
out-of-plane magnetization appears only in a very thin 
region around the vortex core with about the size of an 
exchange length (typically about 10 nm for magnetically 
soft materials'^). The vortex state is degenerated with 
respect to the upward or downward magnetization of the 
vortex core (the vortex polarity p = ±1), hence the vor- 
tex polarity can be considered as a bit of information 
in nonvolatile magnetic vortex random-access memories 
(VRAM). That is why one needs to control the vortex 
polarity switching process in a very fast way. 

The vortex polarity switching phenomena were pre- 
dicted originally for the Heisenberg 2D mag nets.^'"^ The 
interest to this problem was renewed after an experimen- 
tal detection of the vortex core reversal in nanodots by 
an excitation with short bursts of an alternating field,* 
which opened a possibility to use the vortex state dots 
as the VRAM. Moreover, this motivated numerous fun- 
damental studies of the vortex core switching mechanism 
itself. 1 

There are two basic scenarios of the vortex polarity 
switching. In the first, axially-symmetric (or punch- 
through) scenario, the vortex polarity is switched due to 
the direct pumping of axially-symmetric magnon modes. 
Such a switching occurs, e.g., under the infiuence of a DG 
transversal field. ^^^^ In the second, axially-asymmetric 
scenario, the switching occurs due to a nonlinear res- 
onance in the system of certain magnon modes with 
nonlinear coupling. Such a scenario is accompanied 
by the temporary creation and annihilation of vortex- 
antivortex pairs.* The axially-asymmetric switching oc- 
curs, e.g., under the action of different in-plane AG mag- 



netic fields or by a spin polarized current, see Ref. 15 and 
references therein. 

Recently the interest to the axially-symmetric switch- 
ing was renewed: using the micromagnetic simulations 
Wang and Dong^^, Yoo et al.^"^ demonstrated that the 
vortex polarity reversal can be realized under the action 
of an alternating perpendicular magnetic field. In this 
case the resonant pumping of the radial magnon modes 
initiates the switching at much lower field intensities than 
by the DG fields. 

We have very recently predicted the possibility of the 
chaotic dynamics of the vortex polarity under the action 
of the homogeneous transversal AC magnetic field in the 
10 GHz range:!* 

B{t)^e,Bosm{2TTft). (1) 

In order to describe the switching behavior we proposed 
in Ref. 18 the analytical two-parameter cutoff model, 
which gave us a possibility to describe both deterministic 
and chaotic behavior of the vortex polarity. 

The goal of the current work is to study in detail the 
vortex dynamics under the action of a perpendicular AG 
magnetic field: we found a rich vortex polarity dynami- 
cal behavior, including the regular and chaotic regimes of 
magnetization reversal. In order to analyze the compli- 
cated temporal evolution of the vortex polarity we used 
here the discrete reduced core model,^'"^'^^ which allows 
us to describe different regimes of vortex polarity dynam- 
ics, including the resonant behavior, the weakly nonlinear 
regimes, the reversal dynamics, and the chaotic regime. 
The reduced core model is another way to treat the dis- 
cretness effects. As opposed to the cutoff model, the core 
model is simpler, hence it allows to go further in analyt- 
ics. 

The paper is organized as follows: The full-scale micro- 
magnetic simulations are detailed in Sec. II. Our diagram 
of switching events demonstrates regimes of the regular 
reversal (single, multiple and periodic ones), intermittent 
and chaotic regimes. In Sec. Ill we describe the compre- 
hensive vortex core dynamics using a simple collective 
coordinate model, which provides all features of the full- 
scale simulations. We propose a way of a unidirectional 
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FIG. 1: (Color online) Nonlinear resonance curves from 
micromagnetic simulations. Insets: (a) The temporal 

evolution of the polarity, (b) the FFT spectrum of the 
vortex polarity for Bq = 9mT, / — 13.2 GHz during 
15ns. Arrow indicates the pumping frequency. 



switching controlled switching in Sec. IV. In Sec. V we 
state our main conclusions. In Appendix A we derive the 
reduced core mode. We use the method of multiple scales 
to perform a weakly nonlinear analysis of the analytical 
model in Appendix B. 



II. MICROMAGNETIC SIMULATIONS OF 
REGULAR AND CHAOTIC DYNAMICS 

Nowadays the micromagnetic simulations are the in- 
herent tools for the nanomagnetic research.^" Namely us- 
ing the numerical simulations it was shown in Refs. 16 
and 17 that the resonant perpendicular field forces the 
vortex core to reverse. Here we perform full-scale mi- 
cromagnetic simulations to study the complicated vor- 
tex core dynamics in details. We consider a disk-shaped 
nanoparticle (198 nm in diameter and 21 nm in thickness) 
under the action of the vertical oscillating field (1) using 
an OOMMF^^ micromagnetic simulator with integra- 
tion method RK5(4)7FC. The material parameters cor- 
respond to Permalloy (NigiFeig) with exchange constant 
A = 13 p J, saturation magnetization Mg = 860kA/m, 
zero anisotropy coefficient and the Gilbert damping co- 
efficient a = 0.01. The mesh cell was chosen to be 
3 X 3 X 21 nm. For all OOMMF simulations we use 
as initial state the relaxed vortex with the polarity di- 
rected upward and counter-clockwise in-plane magneti- 
zation direction. 

First of all we examined the eigenfrequencies of the 
lower axially-symmetric spin waves by applying a rect- 
angular pulse with the strength of 30 mT during 100 ps 
to the vortex state nanodisk in the same way as in 
Ref. 16. Using the fast Fourier transformation (FFT) 
of the z~component of the total magnetization we iden- 



tified the eigenfrequency of the lowest symmetrical mode 
/m=o = 13.98 GHz. The next nearest peaks in the 
FFT spectrum correspond to 16.75 GHz and 27.93 GHz 
(Fig. 2). 

It is already known^' that the vortex polarity switch- 
ing under the action of the AC field (1) occurs in a wide 
range of field parameters (the field intensities Bq and field 
frequencies /). The minimal field intensity is reached 
at about the resonance frequency /q. In the current 
study we are interested in the long-time vortex dynam- 
ics, which is accompanied by the axially-symmetric po- 
larity reversal mechanism. In all numerical experiments 
we calculated the polarity and the position of the vortex 
as functions of time: The vortex position R{t) is deter- 
mined as cross-section of isosurfaces Mx{R) = and 
My{R) = 0,^^ and the vortex polarity p{t) is determined 
as the average z-magnetization of four neighbor cells to 
Rit). 

To study in details the temporal evolution of the vor- 
tex polarity, we simulated the long-time system dynamics 
during 10 ns with the time step of Ips for a wide range 
of the field parameters (i?o varies from 10 to 110 mT in 
increments of 10 mT, and / changes from 9 to 19 GHz 
in increments of 1 GHz). The results can be summarized 
in the diagram of dynamical regimes, see Fig. 2a. De- 
pending on the field parameters (Sq, /), one can sepa- 
rate several different dynamical regimes: (i) the absence 
of the vortex polarity switching, (ii) the chaotic polarity 
oscillations, (iii) the regular switchings with frequencies 
depending on the field frequency, (iv) the intermittent 
switchings, and (v) the complex vortex-magnon dynam- 
ics, where the vortex escapes from the origin. 

(i) We start from a weak field: the field intensity is 
not strong enough to switch the vortex polarity; this 
regime corresponds to the linear or weakly nonlinear os- 
cillations of the vortex polarity (marked as open boxes □ 
in Fig. 2a) . The weak pumping of the system (field inten- 
sities So < 5mT, see Fig. 1) causes the resonance at the 
frequency /q . The increase of the field intensity leads to 
the nonlinear dynamical regime. However, if the field in- 
tensity is not strong enough, one has a weakly nonlinear 
regime, which corresponds to the nonlinear resonance. 
Apart from the nonlinear resonance behavior, the strong 
pumping causes the vortex polarity instability,^'' it also 
causes the shift of the main peak in the FFT spectrum 
(see Fig. lb) and the beats in the polarity oscillations 
(see Fig. la). 

Let us consider the case when the magnetization re- 
versal occurs. The switching diagram (see Fig. 2a) has 
two well-defined minima. The first one corresponds to 
the resonant excitation of the radially symmetrical mode 
/g. The second minimum near 18 GHz, probably, corre- 
sponds to the dynamics near the higher resonances^^. 

(ii) The red circles • on the switching diagram (see 
Fig. 2a) correspond to the chaotic polarity reversal pro- 
cess. The typical temporal evolution is shown in Fig. 2d. 
To draw a conclusion about chaotic behavior of the vortex 
polarity, or more accurately, to make quantitative mea- 
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FIG. 2: (Color online) (a) Switching diagram: open boxes □ describe the vortex dynamics without switching and 

circles correspond to parameters where the polarity reversal is observed (the red circles with cross indicate 
parameters where the vortices escape from the origin during the first 10 ns, the red circles • represent parameters 
where the autocorrelation function (2) rapidly decays, for the green circles • the switching process is periodical, and 
the red open circles O correspond to an intermittent process), (b) The pseudo-phase diagram for the chaotic vortex 
polarity dynamics (the time step 1 ns, the filled circles mark polarities p = ±1); (c) the same diagram for the regular 
dynamics, (d) The example of the chaotic process, (e) The example of the intermittent process, (f) The example of 
the regular process, (g) The examples of the vortex trajectories in case of the complex vortex dynamics. Black 

points mark the places of the polarity switching. 



sures of chaotic dynamics, we use two standard ways: 
the autocorrelation function for the temporal evolution 
of the vortex polarity and the Fourier distribution of its 
frequency spectra.^'* 

First, we define the autocorrelation function of the vor- 
tex polarity signal 

^(*') = ]^Ep(*'+^>(*j)' * = (2) 

for the discretized time tj — jtQ with the step Iq — 1 ps, 
with the boundary values assumed as zero. It is well 
known^'* from the correlational analysis, when a signal 
is chaotic, information about its past origins is lost, i.e. 
the signal is only correlated with its recent part: the 
autocorrelation function decays very rapidly, C{t) — 
as t — ?► oo.^"* For a periodical signal, the autocorrelation 
function is a periodic too. A typical example is presented 
in Fig. 3: the autocorrelation function C{t) is aperiodic 
and sharply decays for the applied magnetic field 70 mT 
with the frequency 14 GHz, which corresponds to the 
chaotic dynamics. The autocorrelation for the regular 
dynamics demonstrate the oscillations under the action 
of Bo = 110 mT with / = 18 GHz. 



The second way is to calculate the Fourier spectrum 
of a chaotic signal. A typical FFT signal is presented in 
the Fig. 4. It is distinctive for the chaotic regime that 
the continuous spectrum dominates the discrete spikes 
(one can identify in the Fig. 4 only one discrete spike 
at the pumping frequency). The fitting of such a signal 
demonstrates that the spectrum looks like a pink noise. 

(iii) The regular oscillations of the vortex polarity ap- 
pear in the high frequency regime, see the green circles 
f on the switching diagram (see Fig. 2a). We have de- 
tected the periodical motion of the vortex polarity using 
the pumping frequency 18 GHz with the field intensities 
higher than 100 mT. The main peak in the FFT spectrum 
corresponds to 6 GHz, i. e. it occurs at //3 of the pump- 
ing. Other spikes with decaying intensities appear with 
steps of 6 GHz. We compare the autocorrelation func- 
tions for the regular and chaotic oscillations, see Fig. 3. 
In contrast to the chaotic regime, C(t) for periodic oscil- 
lations exhibits a high periodicity with a slowly decaying 
amplitude due to the finite observation time. 

In order to compare the temporal dynamics of the po- 
larity in chaotic and regular regimes, we calculate the 
pseudo-phase trajectories. The method of the pseudo- 
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FIG. 3: (Color online) Autocorrelation function (2) for 
Bo = llOmT and / = 18 GHz (blue dashed line) and 
Bo = 70 mT and / = 14 GHz (red solid line). The first 
process demonstrates a periodical process with 
periodical C(t) and the second one demonstrates a 
chaotic behaviour with rapidly decaying C{t). 



phase space is usually used when only one variable 
[the discretized vortex polarity p{ti) in our case] is 
measured:^^ the pseudo-phase-space plot can be made 
using p{ti) and its future value Pti-^.!, where the absolute 
value of the time step i^+i — ti affects only the shape of 
the trajectory. In the case of chaotic dynamics, one has 
open trajectories in pseudo-phase-space (p(fi),p(ii+i)), 
see the Fig. 2b. In the regular case, pseudo-phase tra- 
jectories are closed, see the Fig. 2c. Both trajectories 
are shown for the first 10 ns of the dynamics: in the first 
case the trajectory every time makes a new loop in a 
different place and in the second case all loops coincide. 
Below in Sec. Ill we construct the phase trajectories for 
the theoretical model of our system (see Figs. 8c, 8b). 

(iv) The red opened circles O on the switching dia- 
gram (see Fig. 2a) correspond to an intermittent pro- 
cess. The typical example of the temporal dynamics in 
such a regime is plotted in the Fig. 2e: the vortex state 
can retain its polarity for a relatively long time of a few 
nanosecons; after that multiple reversal processess occur 
during 50 — 100 ns. Note that in the vicinity of other 
regimes in the switching diagram we observed that the 
vortex polarity, after a few switches, can be 'frozen' for 
the rest of the observation time. For example, two switch- 
ing events occur during the first 1.2 ns (Bo = 30 mT and 
/ = 14 GHz, see Fig. 2c), after that the dynamical polar- 
ity has only weak oscillations. A similar picture occurs 
for Bo = 70 mT and / = 17 GHz, where after three re- 
versals during the first nanosecond the resulting polarity 
remains negative. Since the reversals occur only at the 
beginning, one can conclude that this occurs because the 
field is not switched on smoothly. 

(v) The last regime corresponds to the field param- 
eters [Bq, /), where the vortex escapes from the sys- 
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FIG. 4: (Color online) FFT spectrum of the vortex 

polarity {Bo = 70 mT, / = 14 GHz): solid line 
corresponds to the numerical data, the dashed line is 
the fitting to the pink noise. 



tem origin on a long time scale (see the red circles with 
crosses ® in Fig. 2a). Typically, the vortex starts to move 
during the first 10 ns. The detailed analysis shows that 
the switching scenario differs essentially from the above 
mentioned one: the magnetization reversal is mediated 
by the transient creation and annihilation of a vortex- 
antivortex pair* (for details of the axially-asymmetric 
switching mechanism see Ref. 15 and references therein). 
Two examples of the possible trajectories are shown in 
Fig. 2g: The trajectory gi corresponds to the chaotic mo- 
tion, which is accompanied by numerous reversal events. 
In the regular regime the vortex trajectory has a smooth 
shape (g2). When the vortex stays in the center of the 
sample, polarity switching is accompanied by generation 
of the radially symmetrical modes. After some time of 
observation, a new 4-fold symmetry occurs around the 
vortex, which was mentioned in the Ref. 17 and linked 
with the square mesh symmetry used in the OOMMF. 
When the vortex moves from the center, the switching 
scenario is changed: a pair antivortex-new vortex is cre- 
ated and the antivortex annihilates with the old vortex. 
The further magnon dynamics becomes unpredictable. 

We performed very long-time simulations (up to 30 ns) 
for all parameters from the switching diagram, where 
the vortex does not leave the disk center (see Fig. 2): 
the vortex motion was found for all parameters with 
/ < 17 GHz. For higher frequencies (e.g., for / = 18 GHz 
and Bo = 100 mT) the small oscillations of the vor- 
tex position were observed only for t > 29 ns. In the 
prolonged simulations (iv) the vortex polarity does not 
change its value during the time of observation in agree- 
ment with the conclusion made above that the field is 
sharply switched on. 

Additionally we computed the switching diagram for 
the low-frequency range (see Fig. 5): the vortex polarity 
dynamics was simulated during 5 ns. One can identify 
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FIG. 5: (Color online) Diagram of switching events for 
the low-frequency range. Notations are the same as in 
Fig. 2. The dashed line shows the border of the 
switching region. 
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FIG. 6: (Color online) Schematic of the reduced core 
model: Thick red arrows (numbers 1,4) indicate free 
magnetic moments and thin black ones (numbers 5, 12) 
indicate fixed magnetic moments with = 0. 



from the plot two local minima (4.5 GHz and 6 GHz), 
which correspond to resonances for fractional frequencies 
(|/q and ^/o)- In the range of parameters, presented in 
Fig. 5, there are three regimes: when the field intensity is 
not strong enough, the vortex polarity is not switched, in 
the same way as in the Fig. 2a. The strong field causes 
the vortex polarity reversal, which corresponds to the 
quasi-static regime and the lower fields cause the escape 
of the vortex form the system origin. We checked the idea 
about the quasi-static regime by computing the thresh- 
old value for the static field, which is in our case 611 niT, 
cf. Refs. 10 and 11. 



III. DESCRIPTION OF DIFFERENT 
DYNAMICAL REGIMES 

To gain some insight to the switching mechanism, we 
need a model which allows the magnetization reversal 
process. It is worth reminding that in the continuum 
limit the vortex states with different polarities are sepa- 
rated by an infinite barrier. In the spin lattice the barrier 
becomes finite^^ and the reversal can occur. It is already 
known from our previous paper^^ that the dominating 
contribution to the switching mechanism is caused by 
the exchange interaction inside the vortex core. That is 
why to describe the complicated polarity reversal process 
we use here the discrete reduced core model, which was 
initially introduced by Wysin for the vortex instability 
phenomenon. Later the vortex core model was developed 
to analyze the vortex polarity switching in Heisenberg 
magnets. 

We consider the anisotropic classical Heisenberg disk 
shaped system with thickness and the radius L, as- 
suming that the magnetization of the magnet is uniform 
along the thickness. The energy of such a magnet with 



the account of the interaction with magnetic field reads 

E = XI ["^" ■ "^"+'5 ^ ^ 

-alMsL,^mr,-B{t). 

n 

Here m„ is a normalized magnetic moment on the site 
n of a three-dimensional cubic lattice with the lattice 
constant uq, the vector d connects nearest neighbors, A 
is the exchange constant, the parameter A S (0, 1) is the 
effective anisotropy constant, and Ms is the saturation 
magnetization. 

In the reduced core approach, we suppose that only 
the four magnetic moments of the first coordinate shell 
can vary, forming the vortex core; all the other moments 
are fixed in the sample plane in a vortex configuration 
(see Fig. 6). By symmetry, all four moments are char- 
acterized by the same out-of-plane magnetization /i and 
equal in-plane phase Tp, which is determined as a devia- 
tion from the vortex configuration. We consider the core 
magnetization /i, which has the meaning of the dynamical 
vortex polarity, and the in-plane turning phase as two 
collective variables. The energy of the model, normalized 
by e = SALzX has the form (see Appendix A for details): 

^ = — - — A^l — /i^ cos-0 — ^hshiLOT, (4) 

where we introduced the reduced anisotropy parame- 
ter A = 2/ (a/5 A), the reduced field intensity h — 
aQMsi?o/(2v4A), the reduced field frequency uj = 
2Tr f Ms/ (ej), and the rescaled time r = ejt/Ms. 

The magnetization dynamics in the reduced core model 
can be described by the following equations (see Ap- 
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pendix A for details): 



/i = A^l — sin?/; + -q — /x^) 

— A/i-\/ 1 — y? cos + — M^) sinwr 
A/icos?/; . r^Asin?/; 



(5) 



ft, sinojr 



where the overdot means the derivative with respect to r 
and 77 is a Gilbert damping coefficient. The ground state 
of the model corresponds to 



Mo 



= ±v/l-A2, 



(6) 
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In terms of the core model two opposite values of /ig 
describe vortices with opposite polarities /ig. 

Let us start our analysis with a system without damp- 
ing, rj = 0. Supposing that the turning phase is small 
enough, \tp\ ^ 1, one can easily exclude ip from the con- 
sideration. In this case the Eqs. (5) correspond to the 
effective Lagrangian 

^ ^.M^ _ ^^^^ ^ ^hsinujT, 

AVl-A*^ 2 

This simplification allows us to interpret the complicated 
dynamics as the motion of a particle with variable mass 
^ in the double-well potential ^ (/i) under a periodic 
pumping, see the inset in the Fig. 7. The linear oscilla- 
tions near the equilibrium state correspond to the har- 
monic oscillations of the effective particle in one of the 
wells; the eigenfrequency of such oscillations is 

ujo = v/l-A2. (8) 

Let us describe the nonlinear regime of the dynamics. 
In spite of the small damping in the system, its value can 
be comparable with the pumping intensity. Therefore 
we consider below the full set of the model equations 
Eqs. (5). To analyze the weakly nonlinear regime, we 
use the multiscale perturbation method^^"'^^ . When the 
field intensity is much less than the frequency detuning 
{h <^ \io — ujq\), we can limit ourselves to a three-scale 
expansion in the form 

3 

^ = Mo + 51 £'>n(ro, Ti, T2), T„ = e"T, 

n=l 

3 (9) 
= ^ e"V'„(ro,ri,T2), Lo = LUo+uj±, 

Using such the expansion, one can derive from Eqs. (5) 
the resonance curve u)±(hi) as the solution of the following 



FIG. 7: (Color online) Nonlinear resonance curve: The 
numerical solution of Eqs. (5) (the dashed line) and the 
analytical solution (10) (the solid line) for the following 
parameters: A = 0.9415, h = 0.0002, = 0.01, 
m(0) = 0.337, V(0) = 0. The effective double-well 
potential see (7), is plotted in the inset. 

equation 

V * 2A2 ; (10) 

+ 772(i_a2) (2-A2)'a^ 

with \a\ being the amplitude of oscillations, see Ap- 
pendix B for details. The typical nonlinear resonance 
curve is shown in the Fig. 7, cf. Fig. 1. 

We analyze the strongly nonlinear regimes solving nu- 
merically Eqs. (5) in a wide range of parameters (w,/i), 
see the diagram of switching events in the Fig. 8a. The 
absolute minimum in the diagram corresponds to the 
switching in the range near the resonance frequency wq. 
Other local minima correspond to resonances at dou- 
ble frequency, 2a;o, and sub-harmonics, wo/2 and wo/3. 
Note that all resonance frequencies are shifted in the low- 
frequency direction due to the nonlinear resonance. 

We classified the dynamical regimes by using the 
method of Poincare maps (15 • 10^ points per map). We 
constructed such maps for each pair (cj, K) where the 
switching takes place. One can separate four oscilla- 
tion regimes related to the corresponding regimes in the 
OOMMF simulations (Sec. II) with vortex dynamics in 
the center of the sample: (i) the absence of switching, 
(ii) the chaotic dynamics, (iii) the regular porlarity os- 
cillations between two polarities i^Oi (iv) the switching 
with final oscillations around one of the points (±/io,0) 
in the coordinates (/i, ')/')■ 

(i) The border between the switching region and the 
region, where field or frequency are not enough for jumps 
between (±^OjO), shows a few well-defined resonance 
minima corresponding to resonances at a;o/3, i^o/2, ^^o 
and luiQ. 
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FIG. 8: (Color online) (a) Map of dynamical regimes in the coordinates oj/ojo — h/rj. -q = 0.002, other parameters 
are the same as in Fig. 7. The red region corresponds to the chaotic dynamics. Three scales of grey indicate (in 
order of growing intensity) the switching process with one stable focus on the Poincare diagram and the final 
dynamics near the upward or downward polarity. The light-blue region corresponds to the dynamics with more than 
one stable focus on the Poincare diagram, (b) The phase diagram (projection) for the chaotic dynamics. The black 
circles indicate the equilibrium polarities ±/io- (c) The same as (b) for the regular dynamics, (d). Example of the 
chaotic oscillations. ti = 186 • 10'^. (c) Example of a process whith the final state near — /io- (f) Example of the 
regular dynamics, including the //3 peak in FFT spectrum, (g) The same as (b) for the dynamics with 5 stable 

focuses on the Poincare map. 



(ii) The chaotic dynamics occurs for the low-frequency 
part of the Fig. 8a and in the stretched region between 
resonances at wo and 2wo- An example of the temporal 
evolution is shown in the Fig. 8d. The projection of the 
phase diagram on the (/i, ip) plane (see Fig. 8b) looks 
similarly to the pseudo-phase diagram in the Fig. 2b: 
the projection of trajectory is not closed and representa- 
tion point makes a lot of windings around both ±/j,o- The 
shape of the chaotic Poincare maps dependens on the fre- 
quency. They show the shape of strange attractors, see 
Fig. 9. Their Cantor structure is ill-defined due to using 
low damping. Note, that they are similar to strange at- 
tractors for the Duffing oscillator^^ (nonlinear oscillator 
with quadratic and cubic nonlinearities in the double- 
well potential). However, the reduced core model has a 
more complicated nonlinearity term; using the mechani- 
cal analogy one can speak about the motion of a particle 
with a variable mass (7) in a double-well potential. 

(iii) The main part of the diagram of switching 
events 8a is occupied by the region of regular dynam- 
ics. The most frequently observed Poincare maps for 
this case contain some number of stable focuses. The ob- 
served numbers are 1,12, 15, 16, 18, 21, 24, 30 and 96. 
The most frequently observed ones are 1 (the grey region 



in the Fig. 8a) and 3 (included into the light-blue region 
in the Fig. 8a). Some of the points with a higher num- 
ber of focuses demonstrate a complicated regular dynam- 
ics in phase space, see Fig. 8g. The analogue of quasi- 
rectangular regular polarity oscillations in OOMMF sim- 
ulations is found in the uja/Z region, see phase diagram 
in the Fig. 8c and temporal evolution in the Fig. 8f (com- 
pare with the pseudo-phase diagram shown in the Fig. 2c 
and the temporal evolution in the Fig. 2f). 



(iv) The analogue of the intermittent switching linked 
with perturbation by the applying of the external field are 
shown by two dark-grey color intensities in the Fig. 8a. 
The final dynamics is an oscillations around upward or 
downward polarity (points (/io,0) and (— /io,0) in the 
phase space projection, respectively). An example of 
the temporal evolution is shown in the Fig. 8e. As in 
OOMMF simulations, such oscillations typically occur 
near the border of the switching region. As in case of the 
chaotic dynamics, the resulting polarity is highly depen- 
dent on field, frequency and integration conditions. 
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FIG. 9: (Color online) Evolution of strange attractors with change of pumping frequency, h/i] — 10.5, other 
parameters are the same as in Fig. 8a. The number P means number of points on the corresponding Poincare map. 
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FIG. 10: (Color online) (a) Number of switchings as function of number of oscillations in the wave train (N) and the 

field intensity (i?o)- Symbols correspond to values iV and So, where unidirectional switching is observed, (b) 
Controlled unidirectional switching by sequence of short wave trains (contained N = 6 periods on frequency 14 GHz, 
duration 0.43 ns) with period 3.2 ns. Gray regions show time of applying external magnetic field, (c) The same for 
core model. Separate pulse contains iV = 3 periods (duration 58.9 in arb. u.), interval between pulses 589 arb. u. 



IV. CONTROLLED SWITCHING 

As it is shown by analysis of the diagram of switch- 
ings events (Fig. 2a), the vortex polarity switching un- 
der the action of the perpendicular resonant field pro- 
duces multiple switchings during a short time compara- 
ble with one period of the acting field. Such a situation 
is not unique among different polarity switching meth- 
ods. So, for the axially-asymmetric scenario, sufficient 
strength of gaussian pulse in the sample's plane pro- 
duces more than one sequential vortex-antivortex pair 
creation and annihilations.^^ Using an in-plane rotating 
field with frequency tUr codirectional with the vortex po- 
larity p {uJrP > 0) stabilizes the vortex in the center of 
the sample. But when cUrP < the reversal occurs in 
a specific range of the field intensities and frequencies, 
above which multiple switching was observed. A similar 



picture was reported in Ref. 30 for the current -induced 
switching. 

But for the further application in contrast to multi- 
ple switching a unidirectional vortex polarity reversal, is 
needed. Because the pumping (1) does not select any 
direction of the vortex polarity, the most natural way to 
avoid the multiple reversal consists in limiting the pulse 
duration. We test an influence of a short wave train in 
the form 



B 



BqB^ sin (27r/t) , t G 



0, 



N' 



.0, 



/, 

otherwise. 



(11) 



where g N is the number of periods of the sinusoidal 
magnetic field in the wave train. We investigate the vor- 
tex dynamics under the action of the field (11) on the res- 
onant frequency / = /o- The vortex polarity is observed 
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during a long time lON/f in order to damp magnons. 

The response of the magnetization to the field 
(11) shows the nonlinear dependency on Bq and N, 
see fig. 10a. When field intensity and period numbers 
are small, there are only small oscillations of the magne- 
tization inside the vortex core. When N becomes larger 
than some critical value, a typical system behaviour looks 
like a few switchings, which also occurs when the field is 
already turned off. A unidirectional switching from up- 
ward polarity to downward is observed for Bq = 30 mT 
and N = 6,7 and some higher fields. We check the con- 
trollability of the discussed switching method by apply- 
ing of series of wave trains. The series of wave trains 
is applied to the relaxed vortex. The time interval be- 
tween trains is varied with steps of 0.5//o in different 
series. For Bq — 30 mT and iV = 6 the first switching 
occurs in 627 ps and the nanodot starts to relax (field is 
turned off at the time 429 ps) . The sequence containing 3 
wave trains allows to get controllable unidirectional vor- 
tex polarity reversal with minimal interval between wave 
trains of 3.2 ns, see Fig. 10b. These time intervals corre- 
spond speed of changing state of such memory cell about 
250 MHz. The core model also gives the same qualitative 
result, see Fig. 10c. 



particle starts to oscillate near the bottom of one of the 
wells. When the pumping increases, there appear non- 
linear oscillations of the particle; under a further forcing 
the particle overcomes the barrier, which corresponds to 
the magnetization reversal process. 

The chaotic dynamics of the magnetization is studied 
for domains walls'^ ^ and current-induced phenomena in 
monodomain nanoparticles.^^"'^'* Very recently the exis- 
tence of incommensurate chaotic vortex dynamics in spin 
valves was reported. ^'^ In our case the chaos enters in the 
vortex polarity switching process due to the periodical 
pumping of the system with two equivalent equilibrium 
states as it happens in a Duffing oscillator. 

The periodic pumping does not select the preferable 
vortex polarity direction which causes multiple switch- 
ings under the action of sufficiently high fields and fre- 
quencies. However, accurate fitting of the pulse dura- 
tion and the time interval between sequential pulses al- 
lows to obtain a controlled unidirectional core reversal. 
Thus, the chaotic, regular and controlled vortex polarity 
dynamics could find applications in physical layer data 
encryption'^^''^'' and memory devices.'^'"' 
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V. DISCUSSION 

The axially-symmetric switching of the vortex core un- 
der the action of periodic pumping very recenltly pre- 
dicted in Refs. 16 and 17. Wang and Dong^^ were con- 
cerned with switching events: for the typical nanodot 
size the switching at the resonant frequency occurs dur- 
ing 600 ps. Yoo et al. ^ computed the diagram of switch- 
ing events where they noticed the existence of resonances 
on double and triple harmonics and shown that the ex- 
change energy becomes higher than the threshold value 
for the vortex core reversal. In this work we study the 
long- and short-time vortex dynamics and propose an 
analytical model which describes the phenomena of full- 
scale simulations. 

Our two-parameter cut-off model explains the switch- 
ing phenomenon in terms of the nonlinear resonance in 
a double-well potential. Such a potential arises mainly 
from the exchange interaction: the presence of two wells 
corresponds to the energy degeneracy with respect to the 
direction of the vortex polarity (up or down) ; the energy 
barrier between the wells becomes higher as the discrete- 
ness effects become less important. One has to stress 
that the switching process is forbidden in the continuum 
theory. In a real magnet the magnetization reversal is 
possible due to the discreteness of the lattice. That is 
why to describe the switching analytically we used the 
cut-off model (the cut-off parameter imitates the discret- 
ness effects) . In terms of the cut-off model^® the switch- 
ing can be considered as the motion of an effective me- 
chanical particle with a variable mass in the double-well 
potential. Under the action of periodical pumping the 
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Appendix A: Reduced Core Model 

Taking into account the explicit form of the mag- 
netic field, (1), in terms of the variables m„ = 

(\/l"^-^cos(/)„, A/r^^m2^sin(/)„,TO„j, the energy (3) 

reads^^ 



cos((/)„- <i)n+s) 
+ \mnmn+s\^ - alMsLzBo sin (27r/i) ^ m„. 

n 

(Al) 

We assume that only four magnetic moments in the vor- 
tex core are free, the other ones are fixed in the vortex 
distribution 



(A2) 



with zero out-of-plane component, where x is the polar 
angle, see Fig. 6. We consider the effective vortex polar- 
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ity as the out-of-plane component of the core magnetic 
moments and turning phase 

= fi, V = 0j - (fly, i = 1, 4. (A3) 

According to the Eq. (A2) the energy (Al) reads 

E = -AalMsL^^iB^ sin(27r/t) - AAL^Xy?- 



16 



ALz ^/l — cos -0- 



(A4) 



After the renormahzation the Eq. A4 takes the form (4) , 
where g = E/e,e = iAL^X. 

The magnetization dynamics foUows the Landau- 
Lifshitz-Gilbert equations 



dm„ dS dm 
= m„, X h rjinrtn x 



At 



drrir 



dr ' 



(A5) 



with 7] being a Gilbert damping coefficient and the 
rescaled time r = e^t/AIg. Substituting Eqs. (A3) into 
the Eq. (A5) we obtain the equations for the (ji, i}}): 



(A6) 



d/i dg , ^^dS 

d7 = 9v;-^(^-'^)a^' 

dV' _ dS 7] dS 
dr dji 1 — y? dip ' 

Substituting Eq. (4) into (A6) we get the Eq. (5). 



Appendix B: Weakly nonlinear analysis 

Let us consider the weakly nonlinear case for Eqs. (5). 
Using the series expansion (9), the time derivative reads 



dt 



n=0 



dr ' 



and the equations of motion (5) break into three pairs of 
equations for the different orders in e: 



1- A2 



Difii + L>oM2 = -Vl - A2/iiVi + 



A2 



(B2) 
(B3) 
(B4) 
yl (B5) 



(B6) 



Vl- A2(^2V'l + ^^^^2) - -^i^l 

6 



- (l-A2)7y2/il, 
D27P1 + DitP2 + Dot/js 



_A2 



(B7) 



2A6 



A4 mM2+2^A^iV^i 



+ \/l - A'^7piip2 + ^3 sinwi — ■q24>i- 

The solution of Eqs. (B2) and (B3) reads 

fii{To,Ti,T2) = A{Ti,T2)e"^'>^'' + A*{Ti,T2)e-'^°^'>. 

Following the Floquet theory2^, one needs to omit all sec- 
ular terms. Thus, Eqs. (B4) and (B5) show A{Ti,T2) = 
A{T2) and Eqs. (B6) and (B7) gives the equation for 

Am-. 



2iVl - A2A' + ^772^1- A2 (2 - A^) A 

_2^±^AM*=A2^e-^^^ ^^^^ 
2i 



A2 



(Bl) By solving the Eq. (B8) one obtains the Eq. (10). 
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